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Abstract 


We report on the lunar and solar measurements used to determine the changes in the ra- 
diometric sensitivity of the Sea-viewing Wide Field-of-view Sensor (SeaWiFS). Radiometric sen- 
sitivity is defined as the output from the instrument (or from one of the instrument bands) per unit 
spectral radiance at the instrument’s input aperture. Knowledge of the long-term repeatability of 
the SeaWiFS measurements is crucial to maintaining the quality of the ocean scenes derived from 
measurements by the instrument. For SeaWiFS bands 1 through 6 (412 nm through 670 nm), the 
change in radiometric sensitivity is less than 0.2% for the period from November 1997 through 
November 1998. For band 7 (765 nm), the change is about 1.5%, and for band 8 (865 nm) about 
5%. The rates of change of bands 7 and 8, which were linear with time for the first eight months 
of lunar measurements, are now slowing. The scatter in the data points about the trend lines in 
this analysis is less than 0.3% for all eight SeaWiFS bands. These results are based on monthly 
measurements of the moon. Daily solar measurements using an onboard diffuser show that the 
radiometric sensitivities of the SeaWiFS bands have changed smoothly during the time intervals 
between lunar measurements. Since SeaWiFS measurements have continued past November 
1998, the results presented here are considered as a snapshot of the instrument performance as of 
that date. 
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Introduction 


The Sea- viewing Wide Field-of-view Sensor (SeaWiFS) was launched on August 1, 1997 
aboard the SeaStar spacecraft (now called Orbview-2). The first images of the Earth were taken 
on September 4, 1997, and the first lunar measurements were made on November 14, 1997. On 
September 9, 1 997, measurements of the sun were initiated using the onboard diffuser. Solar 
measurements have continued on a near-daily basis since then. 

SeaWiFS is a second generation ocean color instrument. As such, its mission was de- 
signed, in very large part, on the lessons learned from its predecessor, the Coastal Zone Color 
Scanner (CZCS). Those lessons are discussed below. In addition, SeaWiFS was developed as a 
data buy (Barnes et al., 1994a), with the detailed design of the instrument provided by the manu- 
facturer. However, the performance specifications included a requirement for direct lunar views 
to monitor instrument stability (Barnes et al., 1994a). In addition, the specifications called for 
either an internal light source or a solar diffuser as an onboard monitor of instrument stability. 

The manufacturer of SeaWiFS, the Santa Barbara Research Center (SBRC), chose a solar dif- 
fuser. That decision has a fundamental impact on the long-term stability monitoring program for 
SeaWiFS, and that impact is also discussed below. 

a. CZCS Background 

The Nimbus-7 CZCS was launched in October 1978. It was the first satellite sensor de- 
signed specifically for the estimation of pigment concentrations in the ocean. The mission was 
designed as a proof-of-concept experiment (Hovis et al., 1980), and the second generation of 
ocean color satellite instruments, including SeaWiFS, has been developed using the lessons 
learned from the CZCS experiment (Hooker et al., 1992; McClain et al., 1992). One of the most 
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important lessons was the need for a continuous comprehensive sensor calibration evaluation ac- 
tivity throughout the mission. The processing of the CZCS data set was complicated by the time- 
dependent degradation of the scanner’s radiometric sensitivity, particularly in the visible bands 
(little degradation could be detected in the 670 nm and 750 nm bands). This degradation was 
recognized early in the mission, but quantification of the degradation rate was difficult to assess. 
Although the CZCS had internal lamps, they did not illuminate the entire optical train (Evans and 
Gordon, 1994). Therefore, changes in the characteristics of the optical components at the input 
aperture of the scanner could not be determined from measurements of the calibration lamps by 
the sensor. In addition, it was difficult to separate changes in the sensitivity of the instrument 
from changes in the outputs from those lamps. 

A number of investigators applied vicarious calibration techniques to correct the CZCS 
calibration. Viollier (1982) used a set of simultaneous in situ surface reflectance measurements to 
adjust the prelaunch calibration gain factors to yield a reasonable comparison. Gordon et al. 
(1983) and Mueller (1985) used field observations of the North Atlantic and North Pacific, re- 
spectively, to estimate the time-dependence of the degradation by assuming that the measure- 
ments were representative for those areas. Hovis et al. (1985) used high altitude aircraft under- 
flights to estimate the top-of-the-atmosphere radiances for direct comparisons with the CZCS to- 
tal radiances. The most comprehensive analysis was conducted by Evans and Gordon (1994) who 
assumed that the normalized water-leaving radiances (Lwn) in low pigment open ocean waters 
should match the clear water values of Gordon and Clark (1981). This analysis provided a de- 
tailed time history of the degradation of the visible bands over the entire CZCS mission. How- 
ever, the method assumes that there is no systematic change in the global ocean over the 8-years 
of CZCS operation and does not address changes in the near-infrared bands. It is a reasonable 
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assumption that 520 nm and 550 run are constant for clear water, but may not be the case at 443 
nm where small changes in pigment concentrations can produce significant fluctuations in 
Lwn(443). 

With the exception of Hovis et al. (1985), the vicarious calibration adjustments are de- 
pendent on the of the particular atmospheric correction algorithm applied because the water- 
leaving radiances are small compared to the radiances due to atmospheric scattering. For such a 
calibration to be effective, it is necessary to separate time-dependent changes in the radiometric 
sensitivity of the instrument, including the near-infrared bands, from changes in the atmosphere, 
notably the atmospheric aerosols. As pointed out by Gordon (1987), this requires frequent and 
independent measurements such as images of the moon or views of the sun through a solar dif- 
fuser. As a result, the SeaWiFS mission was designed to accommodate both of these measure- 
ments. 

b. Measurement Background 

SeaWiFS carries no onboard calibration standards. It has a diffuser panel that is used to 
measure the solar irradiance on a daily basis (Barnes and Eplee, 1996). However, the sun is 
viewed by the instrument in a manner different from measurements of the Earth, and the diffuser 
is not used for Earth measurements. It is an extra element used only to view the sun. SeaWiFS 
carries no device, such as a ratioing radiometer (Palmer and Slater, 1991), to measure changes in 
the diffuser’s reflectance. Thus, using solar measurements only, it is not possible to separate 
changes in the reflectance of the diffuser from changes in the radiometric sensitivity of the instru- 
ment. 
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For the SeaWiFS Project, there is one assumption basic to the use of the solar diffuser. 

The change in the reflectance of the diffuser is assumed to be nearly linear for time periods up to a 
few months. Over longer periods, of one year or more, the changes may be approximately expo- 
nential with gradually decreasing changes over time However, this exponential change can be 
treated as a series of linear segments. Experience with diffusers on previous satellite instruments 
(Frederick et al., 1986; Herman et al., 1990) has led to the theory that diffuser degradation on or- 
bit is caused by the coating of the panel with photolyzed organic materials that are outgassed from 
the spacecraft. This accumulation of organic materials is temporally smooth and does not cause 
step functions in the reflectivity of the diffuser. With this assumption of short-term linear change 
in diffuser reflectivity, it is possible to identify sudden changes in instrument sensitivity between 
lunar measurements. 

The SeaWiFS Project does not, as yet, use the moon as an absolute radiometric standard 
for calibration purposes. The moon is used solely as a diffuse reflector whose surface remains 
unchanged (Kieffer, 1997). The SeaWiFS Project cannot, using its resources alone, determine the 
absolute reflectance of the lunar surface nor its absolute radiance. However, lunar observations 
by the US Geologic Survey in Flagstaff, Arizona (Kieffer and Wildey, 1996; Kieffer and Ander- 
son, 1998) are being used to develop a detailed model of the moon that includes such effects as 
libration and phase angle on the reflectance of the lunar surface. The SeaWiFS Project maintains 
an active collaboration with the US Geologic Survey lunar program. 

SeaWiFS views the moon once a month when the moon is about 7° from full phase. The 
selection of this angle is somewhat arbitrary. Because of the inclination of the moon’s orbit to the 
plane of the Earth’s orbit around the sun, there are months where the minimum phase angle for 
the full moon is 1.5°. However, for every month, the minimum phase angle is 6° or less. The 7° 
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phase angle assures the possibility of at least one lunar measurement per month and perhaps two, 
one with the moon approaching full phase and one with the moon leaving. Measurements at 7° 
phase also maximize the illuminated surface of the moon while minimizing the opposition effect, 
the rapid increase in reflected light from the lunar surface as the phase angle approaches zero. 
Operational considerations, such as a conflict of the lunar measurement with a midnight data 
down link, will require the measurements to be moved on occasion to different phase angles. 
This occurred with the lunar measurement in January 1998, where the lunar phase angle for the 
measurement was changed to about 5.5°. There is a lunar phase change of about 0.8° per 
SeaWiFS orbit. By selecting the SeaWiFS orbit closest to 7° phase, the phase angle for each lu- 
nar measurement should be within about 0.5° of the desired angle. 

Instrument Description 

SeaWiFS is an eight band filter radiometer designed to monitor Earth-exiting radiances 
from ocean scenes. It is the only instrument onboard Orb View-2. The sensor’s instantaneous 
field of view (EFOV) is 1 .6 mrad by 1.6 mrad per pixel, with one scan covering 58.3° either side 
of nadir. SeaWiFS can be set to +20°, 0°, or -20° in the direction of flight to minimize the effects 
of ocean glint on the data. Each measurement is digitized to 10 bits, with a typical measurement 
producing about 600 counts with one count of noise. The results of the prelaunch characteriza- 
tion of SeaWiFS are summarized in Barnes et al. (1994a). 

SeaWiFS consists of a scanner, which contains the optics, detectors, preamplifiers, and 
scan mechanisms, and the electronics module, which contains the signal conditioning, command 
and telemetry, and power supply electronics. The SeaWiFS scanner is illustrated in Figure 1 
Light first strikes the primary mirror, an off-axis parabola, and then is reflected from a second 
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surface polarization scrambler and from the half axis mirror before reaching the field stop. The 
half angle mirror removes the rotation of the image from the scan of the telescope. The half angle 
mirror uses alternating sides on successive telescope scans. After the field stop, the light is colli- 
mated by another off-axis paraboloid and directed to the aft optics assembly. Dichroic beam 
splitters in the aft optics divert the light into four focal plane assemblies, each containing two 
spectral bands delineated by narrowband interference filters in close proximity to the detectors. 
Attention in the design of SeaWiFS was given to minimizing the sensitivity of the instrument to 
polarized light. This consideration is the principal reason for splitting the telescope into two sec- 
tions, each rotating at a different speed. This design minimized the incidence angle of light on the 
mirrors. In addition, the use of a polarization scrambler in the fore optics eliminated the need for 
individual scramblers to remove residual polarization at each focal plane assembly. Additional 
details on the design of SeaWiFS are given in Barnes and Holmes (1993). 

For measurements of the sun, SeaWiFS uses a diffuser assembly mounted to the scanner 
(Barnes and Eplee, 1996). The assembly is designed so that the diffuser is illuminated by the sun 
as the spacecraft passes over the South Pole. The diffuser plate is part of the diffuser housing and 
is painted on the inside behind a diffuser cover. The cover is also painted and acts as a second 
diffuser plate. The diffuser cover has a spring-loaded hinge at its bottom and is held in place by a 
solenoid actuator. When the one-time actuator releases the cover, it rotates out of the optical 
path for the diffuser. As of this writing, the diffuser cover continues to act as the instrument’s 
diffuser plate. The plate and cover both have coatings of YB7 1 paint, which provides a durable 
flat-white coating with proven stability, as demonstrated on orbit by the Long Duration Exposure 
Facility (LDEF). As measured by the manufacturer, the paint is spectrally flat over the measure- 
ment wavelength range of SeaWiFS. 
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In order to provide a system-level measurement of the reflectance of the diffuser, the dif- 
fuser housing was illuminated in the laboratory with a source having an angular subtense similar to 
the sun. The illumination source was a 1000 W quartz halogen lamp placed about 305 cm from 
the inlet of the housing. The lamp’s filament subtended an angle of about 1.45° at the surface of 
the diffuser or about three times the apparent diameter of the sun as viewed from the Earth. De- 
termining the reflectance for incident flux normal to the input aperture of the diffuser housing re- 
quired two measurements. For the first measurement, the light from the lamp was measured by 
SeaWiFS from the diffuser. For the second, the instrument was rotated to measure the reflected 
light from a second diffuser. The second diffuser was made of pressed halon and was positioned 
to illuminate the Earth-view aperture of the instrument. The ratio of the two measurements by 
SeaWiFS was used to calculate the diffuser reflectance at normal incidence. For those calcula- 
tions, the reflectance of the pressed halon diffuser was 0.99 /ti. Two dimensional reflectance ta- 
bles, relative to the value at normal incidence, were determined by rotating the instrument while it 
viewed the halogen lamp via the diffuser (Barnes and Eplee, 1996). 

Lunar Measurements 

SeaWiFS operates in a sun synchronous orbit, crossing the equator from north to south at 
local noon. In normal operation, the spacecraft is maintained in a nadir orientation, using pitch- 
axis momentum wheels for attitude control with a spacecraft pitch rate of 360° per orbit. For lu- 
nar measurements, the rotation rate of the momentum wheels is increased, and the spacecraft is 
pitched in the opposite direction at a rate faster than normal operation. The maneuver is started 
past the South Pole passage and is timed such that SeaWiFS will view the moon as the spacecraft 
Earth track passes the sublunar point. At the end of the maneuver, about 28 minutes later, when 
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the spacecraft again points toward the Earth, the pitch rate is returned to normal. During views 
of the moon, the scan direction of SeaWiFS is such that the instrument scans across the lunar 
surface from west to east in celestial coordinates. 

Since the moon appears to be a stationary object during SeaWiFS measurements, the 
number of scan lines in a lunar measurement depends upon the pitch rate of the instrument and the 
apparent size of the moon. The pitch maneuver causes SeaWiFS to over-sample the moon. 

There are approximately 25 scan lines of the moon in the lunar image, whereas the moon has a 
diameter that is equivalent to approximately 7 SeaWiFS samples. With a scan rate of 6 telescope 
rotations per second, the lunar image is collected in about 4 seconds. 

An image of the moon is shown in Figure 2. It is the image for SeaWiFS band 1 for the 
first lunar measurement (November 1997). The image gives the digital counts for each sample 
after the removal of the zero offset. The zero offset comes from a small, constant, internally- 
generated voltage that ensures the digital counts in the data stream are always greater than zero. 
The top of the image (scan line 1) is north, and the left side of the image is west (sample 1). The 
upper left hand sample is designated as scan line 1, sample 1 (1, 1); the lower right hand sample as 
scan line 33, sample 22 (33, 22). The central core of the lunar image includes all samples with 
counts greater than 1% of the maximum. 

In Figure 2, the maximum is 735 counts. The drop off to zero counts at the top and bot- 
tom of the core of the image is approximately the same. There is no such symmetry on the left 
and right hand sides of the lunar image. This effect is due to stray light in the instrument and has 
been seen in laboratory testing of the instrument (Barnes et al., 1995). The moon is a very good 
target to examine the response of SeaWiFS to bright-to-dark and dark-to-bright transitions in the 
scenes it measures. 
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The values in Figure 2 are given as digital counts. This form of the data gives the simplest 
presentation of the measurements. In the SeaWLFS calibration algorithm, however, the digital 
counts are converted to spectral radiances for use in the analysis of the lunar measurements. 

There are factors in the performance of the instrument, such as the temperatures of the focal 
planes and side-to-side differences in the reflectance of the half-angle mirror (Barnes et al., 

1994b), that are part of the counts-to-spectral radiance conversion for SeaWiFS. The use of 
spectral radiances eliminates these instrumental factors from the lunar measurements. 

For the analysis of the lunar measurements, each scene from each band for each measure- 
ment date (such as the scene in Figure 2 for band 1 in November 1997) is represented by the disk- 
integrated spectral radiance. Prelaunch modeling of simulated lunar images (Woodward et al., 
1993), showed that disk-integrated spectral radiances produce better products than those using 
one, or a few, samples from the central image. In the lunar analysis, the summations (disk inte- 
grations) include all of the samples in each 22x33 sample array. They include stray light and other 
instrument-based optical effects. For the summations, the brightest sample in each image ac- 
counts for about 1.5% to 2% of the total. For Figure 2, the brightest sample contains 735 counts, 
so 35 to 40 counts comprise about 0. 1% of the total. The number of counts in the outer border 
of samples in Figure 2 is zero. For the next outer border of samples, the total number of counts is 
one. The use of a 22x33 sample array allows for the inclusion of all parts of the image without an 
excessively large number of samples of deep space. 

a. Normalizing Factors 

Although the surface of the moon remains unchanged over time, the radiance from the 
moon does not. As a result, there are normalizing factors required for the trend analysis. These 
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factors are based, in large part, on the spacecraft positions calculated by the SeaWiFS navigation 
algorithm. The navigation algorithm also provides geolocated Earth coordinates for the meas- 
urements on orbit. It also calculates the location of the instrument above the Earth’s surface at 
the time of the lunar measurements. The SeaStar platform uses a global positioning system (GPS) 
receiver to determine the instrument’s location. The locations of the Earth and the sun relative to 
the moon are derived for SeaWiFS from a calculated ephemeris as a function of the date and time 
of the lunar measurement (Van Flandem and Pulkkinen, 1979). Using these values, five normal- 
izing factors are calculated. 

The first normalizing factor (ki) is the sun-moon distance. Since the sun is an isotropic 
radiator, the reflected irradiance from the moon varies with the inverse square of the sun-moon 
distance. The sun-moon distance, with the moon at full phase, can be calculated as 
Dsm = D se + R (1) 

where D S m is the sun-moon distance in km, Dse is the sun-Earth distance in km, and R is the mean 
radius of the lunar orbit (3.844x10 s km). The SeaWiFS navigation algorithm calculates the actual 
Earth-moon distance, which is substituted for R in Equation (1). When the moon is farther from 
the sun, it is less bright. Thus, the normalizing factor, k>, gives increased values with increased 
sun-Earth distance, with 



( 2 ) 


where kj is normalized to U, the astronomical unit (approximately 1.496x10 s km). 

The second normalizing factor (k 2 ) is the instrument-moon distance. Since SeaWiFS is a 
radiometer with a small, well-defined field of view, there is no inverse square law effect for indi- 
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vidual samples of the lunar surface. However, this analysis uses disk-integrated spectral radi- 
ances, and the integrated image acts as an irradiance source, with 

= (3) 

where Dim is the instrument-moon distance (in km) and Dem is the Earth-moon distance (in km), 
A is the Earth’s equatorial radius (6378 km), and H is the instrument altitude above the Earth 
(705 km). The SeaWiFS navigation algorithm performs a more sophisticated and more exact cal- 
culation of the instrument-moon distance than that in Equation (3). Since the moon fills fewer 
samples when it is farther away from the Earth, the normalizing factor is larger for larger instru- 
ment-moon distances. This factor is normalized to the mean radius of the lunar orbit, using 



The third normalizing factor (ks) is the illuminated portion of the lunar surface as a func- 
tion of the phase angle. This factor is a linear function of the phase of the moon, with the lunar 
surface fully illuminated at 0° phase, half illuminated at 90° phase, and dark at 180° phase. This 
function is given as 

/,((!) = a, +a,ff (5) 

where ao= 1 and ai = -1/180 deg' 1 . 

Factor k 3 is normalized to the fractional area of the moon illuminated at 7° from full phase, using 


09611 

a 0 +a,6 


( 6 ) 


The fourth normalizing factor (Iq) corrects for changes in the brightness of the moon with 
phase angle, which is a function of the change in reflectance of the moon with phase angle. The 
moon has a non-uniform particulate surface, creating large scale regional variations in reflectance. 
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such as variations between lunar mare and highlands. The nonlambertian change in the overall 
reflectance of the lunar surface with phase angle can be approximated by Hapke’s bidirectional 
reflectance equation (Hapke, 1993). Helfenstein and Veverka (1987) used Hapke’s equation, and 
a set of six empirically derived constants, to provide a curve of disk integrated reflectance versus 
phase angle. That curve is shown in Figure 3a; it is given in 1° increments from 0° to 100°. The 
set of coefficients used by Helfenstein and Veverka (1987) are based in large part on previous 
measurements of the lunar albedo (Lane and Irvine, 1973). We use a quadratic fit to provide an 
interpolation between the data points in Figure 3 a. This interpolation scheme is limited to phase 
angles (0‘s) between 4° and 10°, using the function 

f 2 {0) = b o +b l 0 + b 2 0 2 (7) 

where b 0 is 1.287x10**, bi is -6.702xl0' 3 deg 1 , b 2 is 2.163X10* 4 deg* 2 and 0 is the phase angle. 

The quadratic curve agrees with the values from Figure 3a at the 0.1% level. The normalizing 
factor (1 q) is calculated relative to the value at a phase angle of 7°. It is calculated as 

/,(7) _ 0 09238 (8) 

4 /,(») b.+tf+W 

This normalizing factor is shown in Figure 3b. There are indications that the variation in lunar 
reflectance with phase angle has a wavelength dependence. The normalizing factor used here is 
applied over a narrow range of phase angles, and it is anticipated, without complete assurance, 
that the effect of wavelength dependence on this normalizing factor is small. There is also evi- 
dence that the moon is brighter before full phase than after (KiefFer and Anderson, 1998), an ef- 
fect of 0.5% to 1% in the value of Iq. Five of the twelve measurements in the data set presented 
here were made before full phase. 


14 



The fifth normalizing factor is the pitch rate of the instrument during the lunar measure- 
ment. The faster the pitch rate, the smaller the image in the direction of the pitch. Since the 
spacecraft does not have the use of several of its positional sensors during the lunar pitch maneu- 
ver (primarily its horizon sensors), there is increased noise in its internally calculated pitch rate 
during measurements of the moon. As a result, we rely on the number of scan lines in the lunar 
image to determine the pitch rate. To do this, we find the longest vertical section in the image for 
each band. Using this cross section, we determine the points at which the measurement is 1% of 
the maximum value in the section; this is done via interpolation. As a result, the interval between 
the 1 % response points is not limited to an integer number of scan lines. For each lunar meas- 
urement, the intervals for the eight bands are averaged. Since the distances between the 1% re- 
sponse points range from about 24 to about 27 scan lines (Barnes et al., 1998), the results are 
normalized to a value of 25 scan lines, using 



where Lm is the interval between the 1% response points in the longest vertical section of the im- 
age. In addition, the pitch rate normalization process accounts for changes in the Earth-moon 
distance. 

Each of the five normalizing factors used in this procedure is a fraction containing a refer- 
ence constant. The overall normalizing factor for each lunar measurement is the product of the 
individual factors. This multiplicative factor is applied to the summed lunar radiances for each of 
the eight SeaWiFS bands. It can have values as large as 1.07 and as small as 0.89 (Barnes et al., 
1998). For the trend analysis, this result is further normalized to a value of unity for the first lunar 
measurement (November 1997). 
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b. Lunar Libration 


The phase angle is the most important lunar surface parameter for SeaWiFS measurements 
of the moon. The variation of the integrated lunar radiance with phase angle is much stronger 
than any variation with libration angle. For libration changes, the loss of visible lunar surface from 
one side of the moon is balanced by the gain of visible surface from the other. The libration effect 
derives from the difference in reflectance of the gained surface with respect to the surface lost. 
This is expected to be a strong mitigating factor for the libration effect. As with lunar reflectance 
versus phase angle, the effect of libration is expected to be somewhat different for different 
wavelengths. A detailed lunar model is required to account for lunar libration. Preliminary esti- 
mates (Kieffer and Anderson, 1998) show libration to be a 1% to 2% effect for an individual 
SeaWiFS lunar measurement. The complete lunar libration cycle extends for 18 years, and it is 
composed of many subcycles of much shorter duration. For a set of lunar measurements from 
several months to a few years, libration is not expected to have a major effect on the slope of the 
time series. Rather, it is expected to increase the scatter in the data; however, the overall contri- 
bution of libration to the SeaWiFS lunar time series remains unknown to us. 

c. Trends in the Lunar Measurements 

The time series for the SeaWiFS measurements of the moon, covering the lunar year from 
November 1997 to November 1998, are shown in Figure 4. Each of the eight time series in the 
figure is fitted to a straight line to give a first order estimate of its rate of change. The figure also 
shows horizontal lines with values of 0.99, 1.00, and 1.01 as visual references. The data used to 
create Figure 4 are listed in Table 1 . 
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There was no lunar measurement in August 1998. At the conclusion of the July 1998 lu- 
nar maneuver, the satellite’s attitude control system did not reacquire the Earth properly, causing 
the spacecraft to shut down and causing the loss of a few days worth of data. The correction for 
this problem was not fully implemented in time for the August measurement. For the last three 
lunar measurements, starting in September 1998, the reacquisition of the Earth by the spacecraft 
has been uneventful. 

For SeaWiFS bands 1 through 6, the trend lines in Figure 4 have upward slopes around 
0.5% per year. For band 7, the trend line has a downward slope just over 1% per year, and for 
band 8, the annual change is downward by about 5%. The scatter of the data about the trend lines 
is about 0.5%, as listed in Table 2. The scatter for each band is presented in the table as the stan- 
dard deviation of the data points from the trend line. 

There is a distinct pattern to the results in Figure 4. For example, for each band the fifth 
measurement (March 1 998) is highest above its trend line. This is one of five measurements taken 
with the moon before full phase. Since the moon is between 0.5% and 1% brighter before full 
phase than after (Kieffer and Anderson, 1998), this effect adds to the scatter in the measurements. 
The band-to-band similarities in the scatter about the trend lines are an indication that the pattern 
is not an instrumental effect. 

The SeaWiFS geophysical algorithms use the ratios of spectral radiances from the instru- 
ment to derive its ocean data products. These band ratios use relative differences between bands, 
rather than relying on their absolute values. The ratio of band 7 to band 8 is used to derive the 
aerosol radiances in the other bands (Gordon and Wang, 1994), and an accurate relative calibra- 
tion of these bands is essential. The derivation of surface chlorophyll concentration from ocean 
color measurements is simple. The water-leaving radiances from the ocean in the green portion of 
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the spectrum do not change with chlorophyll concentration (Hooker et al., 1992). SeaWiFS band 
5, at 555 nm, measures in the green. However, the water-leaving radiances in the blue bands vary 
inversely with the chlorophyll concentration in the surface waters, since chlorophyll absorbs in the 
blue. The SeaWiFS geophysical algorithms use measurements from two blue bands, band 2 at 
443 nm and band 3 at 490 nm, to provide the blue spectral radiances. From the blue-green color 
ratios, diffuse attenuation and ocean chlorophyll amounts are determined (Mueller and Trees, 
1997; O’Reilly et al., 1998). The use of band ratios also reduces the effects of factors common to 
the measurements from both bands - the effect of surface glitter, for example. 

The use of band ratios can also be applied to the SeaWiFS measurements of the moon. 

For the lunar trend analysis, we have normalized the results for each band for each month by di- 
viding by the average value for bands 1 through 6 for that month. These average values are listed 
in the rightmost column of Table 1. The normalization to the average of bands 1 through 6 re- 
duces the effects of incorrect normalizing factors common to all of the bands - such as an imper- 
fect correction for lunar phase angle. Figure 5 shows the trends in the spectral radiance ratios for 
the lunar measurements relative to the average for bands 1 through 6. As with Figure 4, these 
ratios are normalized to unity for November 1997. 

The trend line for band 5 is very nearly flat, indicating that the normalization used for Fig- 
ure 5 is equivalent to normalizing by band 5, as is done in the SeaWiFS chlorophyll algorithm. In 
Figure 5, the trend lines for bands 1 and 6 are slightly negative, the trend lines for bands 3 and 4 
are slightly positive, and the trend lines for bands 2 and 5 are nearly flat. For these six bands, 
there is no wavelength dependent pattern in the trends. The trends in bands 7 and 8 remain the 
same as those in Figure 4. 
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As shown in Table 2, the scatter about the trend lines for the band ratios in Figure 5 is sig- 
nificantly smaller than that for the individual band measurements in Figure 4. We assume that this 
results from a reduction in the scatter from such effects as lunar libration, to the extent that the 
libration effect is independent of wavelength. For the geometric normalizing factors, such as the 
instrument-moon distance, the contribution to the scatter in the trends should be removed nearly 
completely. In Figure 5, the data point for each band for March 1998 (data point 5) lies almost 
exactly on its trend line. Finally, the trend lines for each band in Figure 5 have values very near 
unity for the first lunar measurement (November 1997). This condition is not found in Figure 4. 

In a previous analysis (Barnes et al., 1998), the trends from the first nine lunar measure- 
ments from SeaWiFS were examined. In that analysis, it was concluded that there was a decrease 
in the radiometric sensitivity of band 6 with a rate of 0.5% annually. With the addition of three 
new data points (September, October, and November 1998), the slope of the band 6 trend line in 
Figure 5 is very close to zero. Trend analyses can, and often do, change with the addition of new 
data points. 

d. Nonlinear Trends 

For SeaWiFS bands 7 and 8, there appears to be a reduction in their rate of change for the 
last three measurements, that is, for the measurements in September, October, and November of 
1998. This suggests that the rate of change for these bands could be approximately exponential 
and could approach zero change with time. These data open the possibility of alternate, non- 
linear functional forms to describe changes in the instrument’s radiometric sensitivity. We have 
found that an exponential curve fit, with a small quadratic correction, 

/?,( t) = e (co+c ' ,+Cl,1) (10) 
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fits the changes in bands 7 and 8 very well. Here, Ri(t) is the relative value of the trend line as a 
function of time, t, in days since the first image on orbit, and the units for the constants in the ex- 
ponent are such that the exponent is dimensionless. For each band, the values for the constants, 
Co, Ci, and C2, are derived from a least squares calculation. The trends for SeaWiFS bands 7 and 8 
are shown in Figure 6, along with the semiexponential trend lines derived using Equation (10). 

The values of the trend lines are close to unity at the time of the first lunar measurement. The 
scatter in the data points about these trend lines is given in Table 2. 

The trend lines in Figure 6 work well over the time interval of the measurements presented 
here. However, SeaWiFS measurements did not stop in November 1998. As a practical matter 
for the production of long-term satellite-based data sets, it is important to provide a means of us- 
ing the existing data to predict the radiometric sensitivity of the bands in the future. Ocean color 
measurements are used by the scientific community in near real time, that is, before new informa- 
tion on the rate of change of the instrument is available. Equation (10) suffers a major drawback 
as a prediction device for future instrument performance. Based on an extrapolation using the 
data at hand, the trend lines in Figure 6 predict an increase in the radiometric sensitivities of bands 
7 and 8 starting in early 1999. Such an improvement is contrary to our understanding of the op- 
eration of the instrument. The increase in sensitivity is allowed by the form of Equation (10). 

There are alternate equations that can be used to describe the change in sensitivity for 
bands 7 and 8. Among them is the exponential function 

R 2 {t)=\-d\\-e- d '‘) (11) 

where R. 2 (t) approaches the value of 1-di over time. Using this functional form, the prediction of 
an increase in radiometric sensitivity is not possible. Because of the absence of a quadratic term 
in the exponential of Equation (11), the trend lines from this equation have greater stiffness than 
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the corresponding lines based on Equation (10). As a result, the fitted curves using Equation (11) 
do not match the changes in slope of the current data for bands 7 and 8 as well as the curves in 
Figure 6. 

Both Equations (10) and (1 1) have another practical drawback to their use with data sets 
that expand over time. With the addition of each new data point to the time series, these equa- 
tions recalculate the trend lines for all of the previous data points. This creates an instability, as it 
were, for the data in the archive. Frequent changes to the data set make the data difficult, if not 
impossible, to use. As a result, the SeaWiFS Project uses a set of piecewise linear trend lines to 
track the changes in the radiometric sensitivities of bands 7 and 8. In this procedure, new seg- 
ments are added to the trend line without changing previous values in the data set. About once a 
year, there is a major reprocessing of the SeaWiFS data set. At these times, it is possible to up- 
date the trend lines from the start of measurements onward. 

Figure 7 shows a set of two-piece linear trend lines for SeaWiFS bands 7 and 8. The first 
segment for each band was calculated using a linear regression plus the first nine data points (No- 
vember 1997 through July 1998). The second was calculated using the last three data points 
(September 1998 through November 1998) for band 7 and the last four data points for band 8. 
The incorporation of the July 1998 data point in the calculation of the second line segment for 
band 7 creates an upward slope for that segment. A single measurement can have a noticeable 
effect on the results from small data sets. Figure 7 was derived with the luxury of an lunar year’ s 
worth of data. In practice, linear segments are updated as new measurements become available. 
The use of several line segments has kept sharp changes in slope, such as those shown in Figure 7, 
out of the SeaWiFS data set. In the first half of 1999, a reprocessing of the SeaWiFS data set is 
planned. At that time, the current set of trend lines will be updated. 
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Diffuser Measurements 


Once in each SeaWiFS orbit when the spacecraft is over the South Pole, the rotation of 
the spacecraft causes the sun to rise and set over the diffuser aperture in the direction of the 
spacecraft’s pitch. Because of the inclination of the orbit of SeaWiFS, the incident solar irradi- 
ance also changes angle over the course of the year in the direction on the diffuser that is perpen- 
dicular to pitch. This angle is called azimuth in the nomenclature of the SeaWiFS diffuser (Barnes 
and Eplee, 1996). The azimuth angles for the SeaWiFS solar measurements are shown in Figure 
8a. They range from about +5° to about —5° from the normal to the plane of the input aperture of 
the diffuser housing. Laboratory measurements of the diffuser’s bidirectional reflectance distribu- 
tion function (BRDF) provide a correction for this seasonal cycle in azimuth. Since laboratory 
measurements were not made for all of the SeaWiFS bands, the BRDF correction for band 8 is 
used here. It is based on Table 29 of Barnes and Eplee, (1996) which gives the laboratory meas- 
urements of the diffuser cover, and of the diffuser, itself. For the correction used here, those 
measurements were fitted to a second order polynomial curve with a value of unity at zero azi- 
muth and a value of about 0.95 at 6° on each side of zero. The effect of the correction using this 
BRDF model is shown in Figure 8b. The correction increases the values of the diffuser measure- 
ments at angles where the azimuth angle is different from the normal to the input aperture of the 
diffuser housing, that is, from zero azimuth. The initial solar measurement with the diffuser was 
normalized to unity in Figure 8b. At day five after the first SeaWiFS image, the correction is 
about 0.6%. For band 8 (865 nm), the effect of the seasonal cycle in the azimuth angle appears to 
be nearly eliminated. 

The solar measurements from the eight SeaWiFS bands are shown in Figure 9. The cor- 
rection for changes in the azimuth angle on the diffuser has been applied to each of them. The 
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values in Figure 9 were selected to cover the time series for the SeaWiFS lunar measurements. 
They have been normalized to unity on day 71 . In this regard, the format of Figure 9 duplicates 
that for Figures 4 and 5. For band 1 (412 nm), there is an apparent dip in the trend line for days 
100 to 200 after the first image. This is the time during which the solar azimuth angles on the 
diffuser are greater than zero. We postulate that this dip is caused by an imperfect BRDF correc- 
tion. The size of the dip decreases for bands with wavelengths closer to that of band 8. Since the 
BRDF correction is based on measurements of band 8 alone, it seems likely that the imperfection 
in the correction increases as a function of the difference in wavelength from band 8. For band 7 
(765 nm), the dip in the trend line is small but noticeable. These changes are seasonal in nature, 
and we assume that they will repeat from year to year. This hypothesis will be tested as SeaWiFS 
continues through its second year of operation. With this additional data, it may be possible to 
derive a correction for the repeating seasonal signature in the diffuser measurements. 

There is also a long-term decrease in the diffuser measurements that is separate from the 
seasonal changes. For bands 1 through 6, there is no corresponding decrease in the lunar meas- 
urements (see Figure 5). The changes in the diffuser measurements for bands 1 through 6 is con- 
sistent with the effects of the build up of a coating on the surface of the SeaWiFS diffuser. Previ- 
ous instruments, such as the Solar Backscatter Ultraviolet (SBUV) radiometer onboard Nimbus-7 
(Cebula et al., 1988; Herman et al., 1990), experienced similar effects in measurements with their 
onboard diffusers. For these instruments, and for SeaWiFS, the long-term changes in the diffuser 
measurements have a wavelength dependence: the changes are greatest in the near ultraviolet and 
the blue, and the changes decrease with increasing wavelength. Such a trend is not seen in the 
diffuser measurements for SeaWiFS bands 7 and 8. For these bands, the long-term changes in the 
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diffuser measurements appear to include the effects of the radiometric sensitivity changes shown 
in the lunar measurements, plus the effects of changes in the diffuser, itself. 

Presently, our imperfect understanding of the seasonal and long-term changes in the dif- 
fuser measurements preclude their use for monitoring the radiometric sensitivity of SeaWiFS in a 
quantifiable manner. The diffuser measurements can be used, however, to check for sudden 
changes in the instrument response between lunar observations. We find no evidence of such 
sudden changes at the level of the short-term repeatability of the diffuser measurements, which is 
about 0. 1%. This is the purpose for which the solar diffuser was incorporated into SeaWiFS. 

The analysis presented here also underscores the importance of a detailed characterization of the 
diffuser’s BRDF before launch. 

Concluding Remarks 

SeaWiFS measurements of the moon show changes in the radiometric sensitivities of 
bands 1 through 6 to be small, about 0.2% or less for the lunar year from November 1997 to No- 
vember 1998. The scatter of the data points about the trend lines is also less than 0.2% for these 
bands. In our analysis, each band has been treated individually, since each has its own interfer- 
ence filter and detector/amplifier system. However, in normalizing the trend data to minimize the 
effects of imperfections in the geometric factors, our analysis may have eliminated a change in the 
instrument sensitivity that is common to all of the SeaWiFS bands. For example, there are optical 
components (including the primary telescope, polarization scrambler, and half angle mirror) that 
may have changing properties effecting all of the bands. There are three factors that lead us to 
believe this is not the case. First, the trends for bands 1 through 6 in Figure 4 are all positive. 
There is no sign of instrument degradation for these bands in this figure. Second, there is no sign 
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of a wavelength dependence in the trends for bands 1 through 6 in Figure 5. For changes result- 
ing from the coating of optical surfaces there is generally a wavelength dependence, with changes 
to a greater degree in the blue and changes to a lesser degree in the red. This appears be the case 
for the SeaWiFS diffuser for bands 1 through 6, as shown in Figure 8. However, in Figure 5, 
there is no sign of such a wavelength dependent effect. And finally, there is no sign of long term 
changes in the radiometric sensitivities of bands 1 through 6, based on comparisons with ground- 
truth measurements during the first year of the SeaWiFS mission (McClain et al., 1998). 

For SeaWiFS bands 7 and 8, the decreases in radiometric sensitivity over the lunar year 
have been about 1.5% and 5% respectively. There are signs that both rates of decrease are slow- 
ing over time. We can fit the radiometric changes in these bands to an exponentially-based non- 
linear function. The scatter of the data points about these trend lines for bands 7 and 8 are less 
than 0.3%. The exponentially-based function derives from the same model that is behind self- 
limiting processes, such as radioactive decay, first order chemical kinetics, and light transmission 
through an opaque medium (the Beer-Lambert law). However, there are practical considerations 
that make the use of this function, or other analytical functions, less than ideal. These include the 
need to predict the future radiometric sensitivities of the bands and the requirement that data 
points added in the future do not change the radiometric sensitivities for measurements currently 
in the data set. 

For these reasons, the trend lines used by the SeaWiFS Project for bands 7 and 8 are a set 
of piece wise linear segments. Segments are added by the SeaWiFS Calibration Team as new lu- 
nar measurements become available. They are based on the collective judgment of the team and 
are used to predict future changes in bands 7 and 8. They show changes with time that are close 
to those in Figure 6. A reprocessing of the SeaWiFS data set is planned for the first half of 1999. 
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At that time there will be a reevaluation of the changes in bands 7 and 8 from the start of the 
SeaWiFS mission. 

We believe the decrease in sensitivities of bands 7 and 8 are caused by changes in their in- 
terference filters and not by changes in their detectors or in the circuits that amplify the outputs 
from their detectors. Electronic checks within SeaWiFS show the detectors and amplifiers to 
function normally, in the same manner as those for bands 1 through 6. 

SeaWiFS carries no onboard device to check for changes in the relative spectral responses 
of bands 7 and 8. A common spectral change is possible, since the two bands show structures in 
their spectral response curves (Barnes et al ., 1994b) that indicate the use of similar dielectric 
components in their construction. However, we currently work with the assumption that the 
changes in the filters are due to wavelength independent decreases in their overall transmission or 
to geometric decreases in their functioning surface area. Such a decrease in active area could re- 
sult from the degradation of the dielectric components at the edges of the filters. To date, the ap- 
plication of this assumption has proved satisfactory in the SeaWiFS atmospheric correction. 
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Table 1 . The data points used to create Figure 4. The rightmost column contains the average 

value for the data points from bands 1 through 6 for each measurement date. This average 


Measurement 

Date 

Days After 
First Image 

Band 1 

Band 2 

Band 3 

Band 4 

Band 5 

Band 6 

Band 7 

Band 8 

MS$M 

14 Nov 1997 

71.27 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

14 Dec 1997 

100.83 

0.9976 

0.9978 

0.9989 

0.9986 

1.0000 

1.0008 

0.9991 

0.9953 

0.9989 

13 Jan 1998 

130.39 

0.9951 

0.9943 

0.9934 

0.9931 

0.9915 

0.9886 

0.9844 

0.9754 

0.9927 

10 Feb 1998 

159.19 

0.9988 

1.0001 

1.0008 

1.0014 

1.0008 

0.9993 

0.9953 

0.9845 

1.0002 

12 Mar 1998 

188.89 

1.0067 

1.0081 

1.0091 

1.0082 

1.0075 

1.0068 

1.0018 

0.9844 

1.0077 

12 Apr 1998 

219.75 

0.9908 

0.9928 

0.9937 

0.9946 

0.9935 

0.9905 

0.9848 

0.9635 

0.9926 

12 May 1998 

249.38 

0.9882 

0.9892 

0.9912 

0.9928 

0.9911 

0.9893 

0.9817 

0.9564 

0.9903 

10 Jun 1998 

278.87 

0.9992 

1.0002 

1.0019 

1.0019 

1.0001 

0.9963 

0.9855 

0.9584 

0.9999 

lOJul 1998 

308.36 

0.9978 

0.9982 

0.9989 

0.9985 

0.9964 

0.9905 

0.9787 

0.9488 

0.9967 

5 Sep 1998 

366.31 

0.9999 

1.0026 

1.0058 

1.0046 

1.0033 

1.0002 

0.9880 

0.9543 

1.0027 

5 Oct 1998 

395.73 

1.0016 

1.0032 

1.0059 

1.0060 

1.0053 

1.0026 

0.9902 

0.9550 

1.0041 

4 Nov 1998 i 

425.84 

1.0005 

1.0035 

1.0055 

1.0058 

1.0048 

1.0017 

0.9882 

0.9535 

1.0036 


Table 2. The standard deviations of the data points from the trend lines. These la values are for 


each panel in Figures 4, 5 and 6, 


* 

Band 

^ 

Center Wavelength 
(nm) 

Standard Deviation 
(Figure 4) 

(%) 

Standard Deviation 
(Figure 5) 

(%) 

Standard Deviation 
(Figure 6) 

(%) 

1 

412 

0.45 

0.13 


2 

443 

0.48 

0.09 


3 

490 

0.49 

0.05 


4 

510 

0.44 

0.06 


5 

555 

0.49 

0.06 


6 

670 

0.58 

0.18 


7 

765 

0.64 

0.27 

0.24 

8 

865 

0.84 

0.49 

0.27 
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direction of flight 



AFT OPTICS 
BENCH 


NADIR 



DIRECTION OF SCAN 
{WEST TO EAST) 


Figure 1. The SeaWiFS scanner assembly. The scanner mounts to the spacecraft using the four 
mounting points at the top of the figure. 
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Figure 2. A lunar scene measured by SeaWiFS. The scene is 22 samples wide by 33 scan lines 

long. The scene is a Mercator projection with the lunar North Pole at the top. The values 
from band 1 are given as digital counts after the subtraction of the zero offset. For the 
analysis of lunar data, these counts are converted to spectral radiances using the SeaWiFS 
radiometric calibration algorithm. 
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Figure 3. Disk integrated reflectance versus lunar phase angle. 

a. The disk integrated lunar reflectance from 0° to 100° phase. 

b. The lunar reflectance normalizing factor. It is calculated using Equation (2). The value at 7° 

phase angle is unity. 
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Figure 4. Changes in radiometric sensitivities of the SeaWiFS bands from the lunar measure- 
ments. The time series are normalized to unity for the first lunar measurement in Novem 
ber 1997. The horizontal lines with values of 0.99, 1.00, and 1.01 are visual references. 
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Figure 5. Changes in the SeaWLFS bands after normalization to the average of bands 1 through 6 
The scatter of the data points about the trend lines is substantially reduced, compared to 
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Figure 6. Changes in SeaWiFS bands 7 and 8. The time series for these bands use the modified 
exponential function from Equation (10). 
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Figure 7. Two piece linear fits to the lunar-based trends for SeaWiFS bands 7 and 8. 






Figure 8. Azimuth angle correction for solar diffuser measurements by SeaWiFS band 8. 

a. The azimuth angle for solar diffuser measurements. 

b. The output from band 8 before and after correction. The symbols give the values before cor- 

rection, the curve gives the values after. The correction increases the values of the dif- 
fuser measurements for azimuth angles different from zero. 
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